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We describe and implement a technique for extracting forces from the relaxation of an overdamped 
thermal system with normal modes. At sufficiently short time intervals, the evolution of a normal 
mode is well described by a one-dimensional Smoluchowski equation with constant drift velocity, 
v, and diffusion coefficient, D. By virtue of fluctuation-dissipation, these transport coefficients 
are simply related to conservative forces, F, acting on the normal mode: F — ksTv/D. This 
relationship implicitly accounts for hydrodynamic interactions, requires no mechanical calibration, 
makes no assumptions about the form of conservative forces, and requires no prior knowledge of 
material properties. We apply this method to measure the electrostatic interactions of polymer 
microspheres suspended in nonpolar microemulsions. 

PACS numbers: 82.70.Dd, 82.70.Uv 



The structure and stability of colloidal dispersions de- 
pend sensitively on the interactions of suspended parti- 
cles. An early triumph of colloid science was the the- 
ory of Derjaguin-Landau-Verwey and Overbeek (DLVO) 
P, Q which juxtaposes short-range Van der Waals forces 
and longer-range electrostatic forces to characterize the 
stability of aqueous colloidal dispersions. Yet the role 
of electrostatic interactions in nonpolar solvents has re- 
mained controversial Electrostatic forces between 
surfaces in nonpolar solvents have recently been reported 
for a variety of surfaces in nonpolar microemulsions 
0, 0, Sj- In certain regimes, measured interactions @ 
are identical to the screened- Coulomb component of the 
standard DLVO theory. In others, a novel counterion- 
only double-layer theory is needed to describe observed 
forces [7] . While the reality of electrostatic interactions in 
nonpolar environments has been established, their origin 
and significance remain mysterious. To that end, robust 
methods for measuring interparticle forces are needed to 
bring out the underlying physics. 

A variety of methods have emerged to directly mea- 
sure colloidal forces. The surface forces apparatus 0, Q 
and atomic force microscope 0, Q measure forces me- 
chanically and have respective force resolutions at the 
nN and pN scales. Alternatively, native thermal fluc- 
tuations can reveal interparticle forces. Such methods 
are well suited to real-space imaging and provide force 
resolutions on the fN scale. For weakly interacting sys- 
tems, the potential can be extracted from the equilib- 
rium distribution [Tot 11 1 by inverting the Boltzmann 
equation U/ksT = — lnP eq . Similarly, liquid-structure 
theory or reverse Monte Carlo methods enable the extrac- 
tion of pair-potentials from the pair-correlation function, 
g(r), of stable semi-dilute dispersions of identical parti- 
cles USE!. This method assumes pair-wise additiv- 
ity of potentials [3] and is difficult to implement with- 
out introducing artifacts from confining surfaces. All of 
these equilibrium methods are limited to interactions of 



less than a few fc^T. An alternative approach, due to 
Crocker and Grier 15], analyzes the dynamics of a sys- 
tem relaxing toward equilibrium. Their method, Marko- 
vian Dynamics Extrapolation (MDE), offers the attrac- 
tive advantage of sampling higher interaction energies by 
driving the system out of equilibrium with an external 
force. MDE elegantly identifies the equilibrium distri- 
bution, P eq , as an eigenvector (with eigenvalue one) of 
the experimentally sampled finite-time propagator. How- 
ever, forces are not extracted from local properties of the 
trajectories. Rather, the propagator must be thoroughly 
sampled over the full range of the interaction - from hard- 
core repulsion at short-range to zero force at long-range. 
Furthermore, systematic effects due to sampling errors 
on the calculation of the eigenvectors are hard to quan- 
tify and artifacts from hydrodynamics are difficult to rule 
out. 

In this Letter, we present a simple method for extract- 
ing conservative forces between isolated pairs of colloidal 
particles from the statistics of their trajectories at short 
time intervals. While our experimental apparatus is a 
straight-forward extension of the blinking optical tweez- 
ers introduced by Crocker and Grier [15| . we propose an 
alternative method of data analysis that measures forces 
locally and implicitly accounts for hydrodynamic cou- 
pling. We apply this method to characterize the elec- 
trostatic interactions of polymer colloids suspended in a 
nonpolar microemulsion. 

We measure the electrostatic interactions of car- 
boxylate modified polystyrene latex particles, radius 
a = 600 nm (Interfacial Dynamics Corp.), suspended 
at vanishingly small volume fraction, < 10~ 6 , 
in a nonpolar microemulsion of AOT (sodium di-2- 
cthylhexylsulfosuccinate) in hexadecane. Samples are 
prepared and stored in a low-humidity glove box. A 
glass chamber, constructed from a standard microscope 
slide, three microscope coverslips (No 1.5), and UV cur- 
ing epoxy (Norland 61), holds the sample for optical mi- 
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pair of optical traps is made with the 532nm output of 
a diode-pumped solid state laser (Coherent Verdi V-5) 
using holographi c op tical tweezers [3, 0] ■ To avoid 
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FIG. 1: Statistics of Particle Trajectories at Short Time In- 
tervals (a) Typical images of trapped and free microspheres 
in 1 mM AOT/hexadecane. (b) A histogram of Ar at various 
At is plotted for an average separation f = 2.04 /mi. (c) Fits 
of the mean displacement and the mean squared-displacement 
as a function of time yield two transport coefficients v and D, 
respectively. 



croscopy and micromanipulation. An inverted optical 
microscope (Nikon TE2000) equipped with an oil immer- 
sion lens (100X, N.A. 1.4) images the suspension in bright 
field. Images are recorded on a high speed digital video 
camera (Photron Fastcam 1024PCI) at a frame rate of 
500Hz. Centroid algorithms [16j| . implemented in MAT- 
LAB, locate particle centers to a resolution of about 10 
nm. 

We extract interparticle forces from the statistical 
properties of the trajectories of isolated pairs of beads. 
We use blinking optical tweezers [15[ to repeatedly trap 
and release particles at a desired separation, as shown in 
Figure QJa). This method allows us to efficiently acquire 
good statistics for the trajectories of spheres in unlikely 
configurations. In our setup, described elsewhere 17j, a 



wall effects [21|, [22J, [23[, particles are trapped at least 
10/im from the walls of the sample chamber. Once the 
optical tweezers have set the height and initial separation 
of the beads, we blink the laser using a chopper (Thor- 
labs MC1000A) at a rate 1/r = 20 Hz, with a duty cycle 
of 1:6. While the laser is off, the particles move freely, 
traveling distances up to about 200 nm. This motion 
is a combination of thermal diffusion and drift induced 
by interparticle forces. As the solvent does not absorb 
a measurable amount of trapping light, we expect that 
local heating by the optical tweezers is not significant. 
Furthermore, if any heat were delivered by the laser, it 
would be dissipated within microseconds by thermal con- 
duction. By studying the motions of the particles only 
when the trap is off, we ensure that our measurements 
are insensitive to the details of the interaction of the trap 
with the particles and to any interactions between parti- 
cles due to light scattering in the traps 24 1 . 

We characterize the stochastic trajectories of free par- 
ticles with time-dependent two-particle probability dis- 
tributions. We reduce the trajectories to a list of statis- 
tically independent events, characterized by initial and 
final particle separations, fj and r/, and a time interval, 
At. Each event is assigned to a spatial bin according 
to its average separation, f = (r^ + r/)/2. We count 
the number of events with a particular value of displace- 
ment, Ar = (rf — r"j), for each separation, f, and time 
interval, At. While each event describes particle dynam- 
ics between a single pair of laser flashes, we concatenate 
events across many flashes to get good statistics. The 
resulting histogram, N(r, Ar, At), represents the time- 
dependent two-particle probability distribution, and is 
well-fit by a Gaussian curve. Histograms for different 
values of At at an average separation f =2.04 [im are 
plotted in FigureQJb). The mean displacement, Ar, and 
the variance of the displacement, o\ r , increase linearly 
with time, as shown in Figure QIc). The slope of the 
mean displacement provides a drift velocity, v. Likewise, 
the slope of the variance provides a diffusion coefficient, 
D. We observe these linear relationships at all particle 
separations. 

The separation dependence of the velocity and diffu- 
sion coefficients are plotted in Figure [2] The relative 
velocity is positive and decreases slowly from a maxi- 
mum at the smallest separation, suggesting a long-range 
repulsive force. In contrast, the diffusion is suppressed 
by viscous forces as particles come near contact. At low 
Reynolds number, force and velocity are related through 
the hydrodynamic mobility tensor, v — b • F. The mobil- 
ity, b, generally depends on the size and separation of the 
spheres and the viscosity of the solvent. If the particles 
and solvent have been characterized in separate experi- 
ments, then the mobility can be calculated from existing 
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FIG. 2: Separation dependence of (a) v(r), the mean velocity 
of separation, and (b) D(r), the relative diffusion coefficient 
in I mM AOT/hexadecane. 



theory 25]. However, typical expressions for hydrody- 
namic mobility are complicated by the presence of surface 
charge and counterions [26| . Without a priori knowledge 
of the zeta potential and screening length, these challeng- 
ing calculations become intractable. We sidestep this ob- 
stacle by exploiting our simultaneous measurement of the 
relative velocity and diffusion coefficients. 

To clarify the relationship between the transport co- 
efficients and forces, let us consider a generalized form 
of Fick's first law for a system of N interacting particles 
123 : 



Ji = 



■D« • [V J P + /3(V J C/)P] 



(1) 



Here, Ji is the probability current of the i particle, Dy 
is the diffusivity tensor, P = P(fi, . . . fjv, t) is the prob- 
ability distribution, (3 = l/fc^T, and U{?\,... r/v) is the 
potential energy. The first term in Equation JI]) is the 
current due to coupled diffusion and the second term cap- 
tures drift due to conservative forces. Enforcing the con- 
servation of probability, we arrive at the Smoluchowski 
equation: 



Vi-D^ • [VjP + p(VjU)P] 



(2) 

The Onsager relations (28| demand that, D,j = D^. 
Therefore, can be diagonalized by a set of normal 
coordinates, x\...xn- If the interaction potential is 
a linear combination of contributions from each mode, 
U = J^i Ui(xi), then the probability distribution can be 
separated, P — \\ i Pi(xi), and 



d t Pi = V i -[T> i -ViPi-v i Pi] ) 



(3) 



where vl = — /3D S ; • VjL/j. Here we have restricted our 
analysis to short time intervals so that the spatial de- 
pendence of the normal modes can be ignored. This 
is valid, provided that: Axiv[ <C Vi, a^ Xi v[ <C fj, 
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FIG. 3: Electrostatic forces between charged colloidal par- 
ticles in a nonpolar solvent. Repulsion between carboxylate 
modifed PS at two concentrations of AOT, as noted. Curves 
are fits to screened- Coulomb interactions (inset) Interaction 
energy, U(r)/k B T, implied by a fit at 1 mM AOT. 



AxiD^ <C Di, and a AxiD'i <C where primed variables 
indicate spatial derivatives. Furthermore, the vector na- 
ture of the normal modes can be ignored provided that 
their displacements are small compared to their magni- 
tude: Axi -C Xi and <J^ Xi -C Xi. This leads to a tractable 
one-dimensional form for the Smoluchowski equation: 

Pi = D % P? - ( Vi - D\)P[. (4) 

For D\ < «i, this is solved by a Gaussian distribution, 



P(Ax t , At) = (2Tro- 2 Ax .)- 1/2 eyip 



(Axi - Aa 



, (5) 

where Axi = ViAt, and o\ x . = 2DiAt. Thus, we ar- 
rive at a convenient expression for the conservative force 
acting on the i th normal mode, 



Fi = -U! = k B T 



D, 



(6) 



This equation, a direct consequence of fluctuation- 
dissipation, relates the force to locally measured trans- 
port coefficients without appealing to particular models 
of hydrodynamic or conservative forces. 

This analysis is readily applied to identical colloidal 
particles, where the separation vector, r, is a normal 
mode of the diffusion tensor and U = U(f). The in- 
terparticle forces can be directly calculated by dividing 
the drift velocities by the diffusivities at each separation 
found in Figure [2j The resulting force profile is shown in 
Figure [3j These spheres show purely repulsive interac- 
tions with a maximum measurable force of about 100 fN 
and a resolution of a few fN. The interparticle forces fall 
off slowly as the particle separation increases, with mea- 
surable repulsions out to about five particle diameters. 
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These results are well fit by a screened-Coulomb form, 

where the Bjerrum length Xb — e 2 /AireeoksT. This fit 
returns a screening length k" 1 = 5.0 ± 0.2 /urn and an 
apparent surface potential, \e(/kBT\ = 3.30 ± 0.04. It 
is important to note that the fitted value of |e£| reflects 
the surface potential as seen from long range. This value 
will be smaller than the actual surface potential for highly 
charged surfaces due to non-linear screening near the par- 
ticle surface. 

Interparticle forces vary with the concentration of sur- 
factant. When the concentration of AOT is increased to 
10 mM, the range and scale of the interparticle forces 
drop, as shown in Figure [3] In particular, the screening 
length is lowered to 0.6±0.I /mi and the apparent surface 
potential is significantly reduced |e£/fcBT| = 1.8 ± 0.1. 
Previous work on a related system [6j reported similar 
values of k and £. 

Interaction potentials, U(r), can be calculated from 
these parameters, as demonstrated in the inset of Figure 
[3l The potential is soft and long-ranged, decaying from 
about 100 ksT at contact to < fc^T at about r = 10 
jtxm. We use fitted values of k" 1 and \e(/kBT\ because 
we have found direct integration of the force curve to be 
highly unreliable; by simply varying the size of the spa- 
tial bins, conventional potentials with scrccncd-Coulomb 
forms can be transformed into anomalous potentials with 
long-range attractions. 

We present a method for extracting the conservative 
forces between colloidal particles from the statistics of 
their trajectories. This method requires no separate 
measurements of solvent and particle properties. The 
only supporting measurements are the spatial and tem- 
poral calibrations of the imaging system and the tem- 
perature of the sample. Additionally, our measurement 
is completely independent of specific models for hydro- 
dynamic and electrostatic interactions between particles. 
Our method of data analysis, rooted in the general prin- 
ciples of noncquilibrium statistical mechanics, may be ex- 
tended to probe generalized forces acting on fluctuating 
normal modes of any thermal system, from more complex 
colloidal systems to the internal dynamics of molecules. 

We thank Ian Morrison and Todd Squires for helpful 
discussions and Cabot Corporation for support. 
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